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The Thomas-Fermi QuEirk Model: Non-Relativistic Aspects 
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Department of Physics, Baylor University, Waco, TX 76798-7316 

Non-relativistic aspects of the Thomas-Fermi statistical quark model are developed. A 
review is given and our modified approach to spin in the model is explained. Our results are 
limited so far to two inequivalent simultaneous wave functions which can apply to multiple 
degenerate flavors. An explicit spin interaction is introduced, which requires the introduction 
of a generalized spin "flavor" . Although the model is designed to be most reliable for many- 
quark states, we find surprisingly that it may be used to fit the low energy spectrum of octet 
and decouplet baryons. The low energy fit allows us to investigate the six-quark doubly 
strange if-dibaryon state, possible 6 quark nucleon-nucleon resonances and flavor symmetric 
strange states of higher quark content. 
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I. INTRODUCTION 

The theoretical aim of our Thomas- Fermi ("TF") statistical quark model, introduced in Ref.[l], 
was to explore many-quark baryon states including strage matter pi [H]. It is a useful new tool for 
the quark physicist to quickly assess the possibility of bound and resonance states in preparation for 
much more detailed and expensive lattice QCD calculations. The TF semi-classical model combines 
Coulombic quark interactions with a bag model type spherical confinement assumption. The model 
was seen to be extremely versatile in that a number of widely varying physics scenarios, including 
non-relativistic, extreme relativistic, massive gluon, or color-flavor locking, can be addressed. The 
development in Ref.jT], however, was limited to an investigation of a single quark wave function 
representing an equal number of A'^^ mass-degenerate flavors. We would like to build upon our 
previous results to show how wave functions for unequal numbers of degenerate flavors or non- 
degenerate masses can be determined though the coupled differential equations. Although the 
model is designed to be most reliable for many-quark states, we find surprisingly that it may be 
used to fit the entire low energy spectrum of baryons. In order to do so, the model had to be 
extended to include an explicit spin splitting term. We introduce spin as a generalized "flavor" 
in our non-relativistic model. This change keeps the quark differential equation unmodified and 
just induces an overall scale change in the energy expressions. The non-relativistic nature of the 
model is a shortcoming, but we view our work as the appropriate initial step toward developing a 
more realistic relativistic version of the model. However, we note that the quark mass parameters 
dominate the overall baryon masses and that the overall fits are extremely good, so the model is 
self-consistent. 

In the next section we will describe some of the changes in the formalism necessary to intro- 
duce quark spin as a generalized flavor. We will then present the differential equations and the 



generalized expressions for the quark energies in Section III for states with two different wave 



functions. More of the formalism for the consideration of spin will then be presented in Section 



IV We will describe the computer program that solves the TF wave function differential equations 
and present our fit to the baryon spectrum in Section [Vj Note that this model is designed only to 
look for true many-quark bound states and resonances rather than bound nucleonic states like the 
deuteron. In this sense, the model complements calculations of loosely bound states as in Ref.[5]. 
As an application, we will attempt such a search in Section VI A| for the hypothesized if-dibaryon 



state [1], make connection to nucleon-nucleon resonances in SectiorVIB, and apply the formalism 



to high multi-quark strange states in Sectior VI C 



3 



II. MATHEMATICAL PRELIMINARIES 

Let us recall how the TF quark model is formulated. The original flavor and color number, 
nlipp), and non-interacting quark kinetic energy, £-{pf), densities were 

and 

where {pf)1 is the Fermi momenta and where the "/" superscript stands for flavor and the "i" 
subscript stands for color. Notice the explicit degenerate spin factors of 2. Such an assumption of 
degeneracy is excellent for atomic physics, marginal for nuclear physics, but unrealistic for hadronic 
physics. The addition of a spin term in the TF Hamiltonian requires the introduction of non-equal 
spin densities. Instead, we will define new quantities distinguished by caligraphic superscripts, 

' ~ J {2T:hf 2m 2Q^mm ' 

where the spin specification is now included in an augmented fiavor-spin index, X. After variation 
of the densities in the energy functional, Eq.(21) of Ref.[l] still holds: 

Ml! = ^ 3x|g' i Nx-l r-- 3 Nx{r') ^ ,3 . n^{r') 

2m^ (3A-l)l Nx J \r-r'\ f^J \f-r'\ 

However, the previous relationship between the TF spatial wave function and the density. 



ra 

is now 



fir) = — ^(3vr2n^(r))2/l (5) 



fir) ^ -^ien'Nfr))^'. (6) 
2 X 

Rather than change the TF differential equations it turns out that all one need do to accommodate 
the new approach is to make a change in the overall spatial scale factor. The previous definition 

2/3 



Rx, R 



2 X 



3ttA 



(7) 
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is now replaced by 



r = TZx, IZ 



2 X %as 



2/3 



equivalent to a simple change in the underlying scale factor: 

a 22/3a. 

The form of the differential TF equations in Ref.[l] and the wave function normalization, 

Nx 



(8) 



(9) 



(10) 



are unchanged from before. 



III. GENERALIZED MODEL FOR TWO WAVE FUNCTIONS 



In Ref. \X] , a non-relativistic model with equal numbers of Nf degenerate flavors was considered 
in a case study. However, in order to do realistic spectrum calculations, we have to generalize to 
unequal numbers of flavors. In this section, we explain how to solve for systems with 2 separate 
masses or systems with two equal masses, but different flavor amounts (numbers). We will see this is 
equivalent to the introduction of two quark flavor-degeneracy factors, (-/V/)i = gi and {Nf)2 = g2- 
The differential form of the TF quark equations for this case, which involves two distinct wave 
functions, is 



A 



ax- 



1 Nx-l 



(/^(^))' + E (/^(^))' 



Here ax is the mass ratio of the lightest mass quark, mi, to quark flavor X 

ax 



mx 



(11) 



(12) 



and Nx is the number of quarks with generalized flavor X. (Note we are using subscripts on A'j, 
OLX and mx so that they will not be mistaken for powers. This is in contrast to Ref.jl] all these 
quantities were denoted with superscripts.) For g\ flavors with equal number A'^i, considered the 
state with the larger radius, x\ > X2, and 52 flavors with equal numbers we have 



giNi + g2N2 = 3 A, 

for quark quantum numbers, and the generalized two wave function equations are: 



ai 



a2- 



dx"^ 
dx"^ 



A 


1 


ViVi-i 


A- 


1 y/x 




A 


1 


'(N2-1 


A- 




[\ N2 



+ 92-1) {f\x)f'+g, ifHx)f' 



(13) 

(14) 
(15) 
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We assume there is a linear relation between f^{x) and /^(x): f^{x) = kf'^{x). The equations (14) 



and (15) become: 



aik- 



a2- 



dx'^ 

d 'fjx) 
dx"^ 



A 1 

A 1 

~A^Vx 



k'/'+92] {fix)) 



3/2 



N2 - 1 
No 



+ 92 -A 



{fix)) 



3/2 



(16) 
(17) 



To make these two equation consistent with each other implies that 



ai5i 



a2 



Ni-1 



+ 51 



1^ A;3/2 



+ ai 



N2 - 1 

iV2 



+ 92-1] k- Q252 = 0. 



We solve this condition numerically for the k value for each particle parameter set. 
Remember we have two normalization conditions: 

J^^'^{fix)f'dx = iN2/3A), 
J^'^{fix)f'dx = iN,/3A). 



(18) 



(19) 



With Equation (19), we can express our normalization conditions in the form of boundary condi- 
tions: 



dx 



fix) 



\X2 



91 



N2-1 I „ 1 
+ 52-1 



N2 

3Aa2 '■ 



(20) 



The fix) = kf'^ix) condition must be consistent with the normalization conditions. Because 
of the attractive nature of the system, it is a natural assumption that the wave functions are 
compacted around the origin. That is, we assume no voids in the particle interiors. Although this 
is a natural assumption, we have not proven these configurations have the lowest energies. 

With these modified boundary conditions, we can derive the expression for kinetic, T, and 
potential, U , energies: 



E 



12 (3^A\^''^ |c/^ X |a. 



Xj 



{f{x)) 



5/2 



-dx 



12 /Stt^^ 



57r 



2 J 



51"! 



{fix)) 



5/2 



-dx + 5202 



-2 {fix)) 



5/2 



-dx 



x 



(21) 
(22) 



12 fSTTAy/' |ff^x|a , 
5tt 



2 

A- 1 



5^ d/2(x)| N2 
3A 
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dx 



a2 



kgiai 



5 '-j^9i<^i dfix )^ ^ 4 a\5/2 , 4 ^^^^ ^^^5/2 

'7^+^^i-l dx 



1 + ^5202 (/^(a:^2)) \/x2 + ^51"! (/^(a^i)) ^/xT 



(23) 



6 



U 



a 



- 1 



^4/3 



(fix)) 

ax- ^ 



dx'Vx'fix') + 



dx\/xf^ (a 



•max J {f\x')) 

\fx' 



4/3 



2 /3^^\^/^ X |as 



vr 



dx'y/^'f'\x')+ I dx^/xf{x) J dx' ^ 

12 gig2kai \ 



4/3- 



(24) 



N2 df{x) 
3 A dx '"^^ 



12 

y 



515201 N2 



A-1/3 ,r 

4 Q'lgi iV2 ,3/2 



giai 



AT r A-1/3 

N2 5 —4^0151 



3A 7^+g,-l 



Ni-l 



, 4 . 0/ ^^5/2 / 4 ,1 x5/2 ^ 

+ ^5202 (j (3^2) j V2;2 + ^fflOl (/ (Xljj y/Xl 

Of course, one is not limited to two different wave functions in this model, and the above 
equations can be generahzed. Such a generalization is not necessary for the low energy baryon 
spectrum. 



5'i52ai N2 
1 3^ 



(25) 



IV. MORE ON SPIN 



As outlined above, the traditional TF model simply assumes a spin degeneracy factor of 2 in 
Eqs.Q and ([2]). We initially attempted to develop such a model based upon Ref.[l]. In this 
traditional approach there is an intrinsic splitting between states like the nucleon and (or 
A^) already just from the different TF flavor wave functions, exclusive of spin, with the correct 
ordering of states. That is, there is a degeneracy splitting from the different quark flavor content, 
uud versus uuu, for example. However, the nucleon-delta splitting turns out to be much smaller 
than the actual splitting for reasonable parameter values. In our attempts to fit the low energy 
baryon spectrum, we found that our and mi parameters (the latter sets the overall scale) were 
driven to unrealistic values to try to account for such splittings. The model was incapable of 
producing a realistic low energy spectrum. It is necessary to add an explicit spin interaction term. 

There is a problem of course because products of the spin "up" and "down" (usually taken along 
the z-axis) in the semi-classical TF model do not combine to form "good" total angular momentum 
states. Non-relativistic baryon states such as the proton and neutron have total angular momentum 
J = 2 and magnetic quantum number, mj. To explain our approach to incorporating spin, consider 
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the appropriately symmetrized proton flavor-spin ^ = 5 j J = 2 ^^^e function: 

\P, + >= \uud > (2| + +- > -| H h > -| - ++ >)/(3V2) + cyclic permutations. (26) 

The TF quark model is incapable of reproducing this linear combination. Instead given this wave 
function the TF model simply asks the probabilities are of certain configurations. In the proton 
world, the possible spin up configurations are: 

The proton is then in the configuration: 

o o 

(We assign no meaning to the flavor or spin sequential ordering.) It is the configurations which have 
distinct masses in the TF model. Thus, the TF spin model deals with probabilities of configurations 
rather than spin amplitudes. We assume the mass of the physical state is the probability weighted 
average of the configurations. Using the flavor degeneracy factors introduced in the last section, 
these two conflgurations are classified: 

{uud) : gi = l,Ni = 2; 52 = 1, -/V2 = 1- 
{uud) : gi=3,Ni = 1; N2 = 0. 

That is, the first configuration has two identical particles (the u quarks) and a second generalized 
flavor (d), whereas the second conflguration simply has 3 non-identical but mass-degenerate quarks, 
with no second set of particle labelings necessary; note how the different spin u quarks are treated 
as different flavors. We will have more to say on the violation of rotational symmetry later in this 

section. 

The classical interaction of a magnetically charged particle in an external magnetic field is given 

by 

{<''nij = -m-Bj, (27) 

where fhi is the magnetic moment of particle i and Bj is the external magnetic field of particle j. 
A magnetic dipole field is given by 

3f(r-m) — m Svr ^ ^, . 
B'^ = ^-/ + —rn5{r), (28) 
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where 5{r) is a Dirac delta function. In our case, the spherical symmetry of all integrations require 
the interaction to be of the form, 

Svr 



class\ 
m Jij 



-fhi ■ fhj S{f). 



(29) 



These considerations lead us to postulate the form of the color magnetic interaction in our model: 

Stt 



(30) 



where the superscripts I and J- are flavor labels and i,j are particle number labels which will take 
on values from 1 to Nx, Nj. We define 



7i 



2m,- c 



(31) 



where g is the gyromagnetic factor, which takes on a value of 1 classically and a value of 2 to lowest 
order in QED perturbation theory for the electron and muon. Note that the specialization to spins 
along the z-axis is appropriate for our individual spin basis. For two quarks. 



Qi ■ Qj 



1 5^, same color 

— I^f^, different color 



(32) 



The treatment of potential energy terms is explained in Ref. [Ij. We sum over same flavor and 
different flavor contributions, using a single particle normalization for the particle densities {n^{r) = 



3nP'{r)/Nx). We then average over the color interactions in Eq.(30) using (32) and the color 
probabilities, Pij. The result for the magnetic spin interaction energy, Em, may be written: 



En 



Svr 9 X 

Y {3A - 1) 



2 f Nx-1 
2Nx 



Nx Nj 



T ^ 'T 1 — 1 



NxNj 



(33) 

/ \ / V-T / V T / 

X<J i=l j=l 

Switching to the TF wave function, f^{r) defined in Eqj6j and using the x variable defined in Eqjs] 
now gives 



Er, 



IQhc 



37r2(3^-l) a 

Ni Nj 



dx 



{f{x)f 



X 



T ^ T A 1 ^ 1 *i 



dx 



1 (/^(^))'/' (/^(x))3/2 



X Nx 



Nj 



X<J i=l j=l 

where the unit normalized normalized spins have {S^)j = 1 for spin up and {S^)j 
down along z. 



(34) 
T for spin 
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As pointed out above, the introduction of spin labeling as an extended flavor attribute has the 
obvious shortcoming of loss of rotational symmetry. One way of stating the issue is that there 
is a conflict between the configurations used to determine flavor content, which use an individual 
spin basis with good values of J|, where i labels the particles, and the rotationally invariant total 
spin states which have good total angular momentum quantum numbers. This problem will not 
affect multi-quark total spin or 1/2 states. For example, the proton m = zt^ states both project 
unambiguously to the configuration ^uud + ^uud. However, the spin | m = | state, 

^(1 - ++ > +1 + -+ > +1 + +- >), 

produces a different mass than the m = | state when projected into the individual spin basis 
states. The m = \ state corresponds to a (71 = l,A^i = 2; 52 = 1,^2 = 1 configuration in our 
model, whereas the m = | state corresponds to g = \^ N = 'i. These have different masses. 

We believe the best projection is always performed in a maximum m states, | J, m = J >. The 
reason has to do with what we will call the "maximum compatibility" of the two sets of basis states. 
As an example, let us consider the spin 1 1, | > 3 quark states. Here there is only one way to build 
the state, namely || >i || >2 >3 in the product spin basis. This state is also an eigenstate of 
total and Jz- The expectation value of the commutator of with the individual Jz vanishes: 

<\,\\[j\j'/^']\\.\>=^- (35) 

This commutator does not vanish in the J = |, m = ^ state. As a more general example, consider 
the J = 0, 1, 2, 3 states produced by coupling 6 quarks. The J = states will have a rotationally 
invariant projection, whereas the J = 3, m = 3 state will be maximally compatible with the 
states with good individual jI''"'^ since the commutator will again vanish. The maximal m states, 
J = 2, m = 2 and J = 1, m = 1 will also be the most compatible states to use with J = 2 and 
J = 1, although here the commutators will no longer vanish. 

V. GROUND STATE BARYON FITS 

In order for the model we are presenting to be predictive, we need to fix the phenomenological 
parameters. These 5 parameters are: 



10 



B: "Bag" constant 
mi: light quark mass 
nis: strange quark mass 
Og'- strong coupling constant 
g: color gyromagnetic factor 



First, we must understand that the model itself is not well designed to fit low quark number 
states, just as the atomic version would be poorly constructed to fit the low atomic binding energy 
states. In fact, applying this model to a "gas" of three quarks would seem to be impossible. 
However, just like the atomic model, an unreasonable aspect of the mathematics of the TF model 
is that such a fit seems to be entirely reasonable. An aspect that helps is the fact that in the 
hadronic case there are both the octet and decuplet ground states energies available which may 
be used to help make the fit more robust; that is, there are more states available than parameters 
which need to be fit. 

There are 13 types of base configurations and associated TF wave functions that are necessary 
to fit all the octet and decuplet states. Table 1 lists the wave functions present in the ground state 
hadrons as well as the particles which partly or wholly share that wave function. For four of these 
wave functions there are both light quark (u and d) and strange quark combinations present. 

In all cases listed in Table 1 except the ones with either 5(1 = 3, A'^i = 3 or A'^s = 3, the 



consistency condition, Eq.(18), does not formally have a solution when the g value or N value is 



equal to one. These cases are indicated in the Table with the designation l"*" to indicate that the 



solution to ( 18 ) is defined by an infinitesimal approach to that parameter from above. For example, 
in the case of the /2 wave functions, which contribute to the energies of the P, and A"*" states, 
we have the robust numerical condition, 

= 1.23162/^(x), (36) 

for the relationship between (/lO ^^'^ {fk) < x < X2 in the limit A^2 — > 1 from above. 
We use this limit to define the value of the model for A^2 = 1 exactly. (Note that f'^\x) = for 

X2 < X < Xl.) 

Table 2 gives both the base particle configuration, to be used in conjunction with Table 1, as 



well as the actual formulas for the spin splittings from Eq.(34) above. Besides the proton and 
neutron, note that a number of particles share the same base particle energies before magnetic 
spin splitting energies are added, such as the A++ the A+'°, the S+ the S*+ the 5° ~, the 
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'-, as well as the A, E° and We list only the A++, A+, E+, E*+, 5°, and H*° particles in 
the Table because of these degeneracies, although the A, Tj^ and S**^ all get different masses after 
spin interaction. 

There is a subtlety concerning the H*^'" and H*'^'^ particles which prevents us from including 
them in our numerical results. All these particles contain two strange quarks and a light quark. 
In our numerical simulations we find that the two strange quark wave function, /| or /|, has a 
smaller radius than the single light quark, /g or fj, which in Table 1 is characterized as either 



= f^, A^^i = 1 or = 1, A^i = 1^. This means that the right hand side of Eq.(14) is zero for the 



region X2 < x < xi, where f'^{x) = 0, and that the light quark is "free" with no restoring force. 



The equations can still be formally solved and the normalization condition for in Eq.(19) can 
still be fulfilled, but we reject this solution as unrealistic. Thus, although the H*^'" and H**^'" wave 
functions are formally calculable, we don't believe the model can not give correct wave functions 
in this case, and they are omitted from our numerical results. 

A computer program in Mathematica has been developed and the parameters which fit the low 
energy baryon spectrum have been determined. Using a grid search algorithm we have found a 
prehminary best fit with the values: ^^"^=80 MeV, mi=300 MeV/c^, = 520 MeV, = 0.44, 
g = 7.31. These values are phenomenologically appropriate and reasonable. The bag constant 
is slightly lower than in other studies, which means the non-relativistic fit is actually more self- 
consistent. The strong coupling constant is consistent with a bag model designed for heavy-light 
systems [6j. The gyromagnetic factor, g, is large compared to the electrodynamic case (~ 2), but 
not unreasonably so. We do a global fit over the 9 particles listed below. Our results so far for our 
mass fit in MeV/c^ units is (experimental numbers in parentheses): 



Spin^: 

P,N: 966 (939) 
A: 1135 (1116) 
E+-: 1214 (1189,1197) 
SO; 1169 (1192) 



12 



Spin |: 

A++'-: 1221 (1232) 
A+'O; 1205 (1232) 
E*+-: 1416 (1385) 
1391 (1385) 
O": 1702 (1672) 

One can see that the model overestimates the strong isospin breaking effects from the particle 
wave functions for the differently charged A, S and S* particles. We note that the A"*""*"'" and A^''^ 
particles actually come from different base configurations; similarly for the and E° as well as 
the S*+'~ and 11*°. There should also be such isospin breaking effects in nature from the exclusion 
principle, but clearly they are not as large as seen here. Although the fit is not as good as for a 
typical non-relativistic quark model, we would not expect it to be. As explained above, the known 
ground state baryons badly violate the many quark semi-classical assumption. Nevertheless, we 
are encouraged by the overall reasonableness of the fit. 

Table 3 gives a full accounting of the various particle energies in our model. Notice that the final 
particle energies involved are consistent with the non-relativistic assumption partly there are large 
cancellations between the kinetic and potential energies. Also notice that the radii are relatively 
large compared, for example, to measured charged particle electromagnetic radii. This is a result 
of the fit with a smaller value of B^^^ than most standard phenomenologies. This is also partly 
responsible for the largely non-relativistic nature of the fits. 

We now turn to try applying the phenomenology developed here to quark configurations in a 
number of sectors. 
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Base configuration 


Particle(s) 


Wave function 






!i-2 




<■>- 






P,N 


/i 
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-- 


-- 


-- 


- 




P,N 
A+ 




1 


2 


1 


1 ' 


-- 


-- 






h 


1 


3 


-- 


-- 


-- 


- 






fifi 


2 


1 


-- 


-- 


1+ 


1 




A 


















s+ 


fiifi 


1 


2 


-- 


-- 


1 


1 ' 






fit fa 


1+ 


1 


-- 


-- 


2 


1 


^ + J ^ ^ ^ 

-V -V ' ill*, ' .s ' u' 
























1 
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1 


2 
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h 










1 


3 



Table 1: Configuration wave functions. 
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Base particle configeration 


Magnetic Energy splitting C ^ ^^^^^ -i)^ {f ) 
















sr-r}..,(-/,,/';'')..,,(/.'f^./.«'f'")} 




















-~"--{(/-'f>-(/-''''i''"')} 




^•■=.p.(/.<f);.(/.(/'f^)} 







Table 2: Spin energy contributions. 
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Particle 


Raclius(f) 


Rest 
mass(MeV) 


T(MeV) 


U(MeV) 


Spin(MeV) 


Volume 
Energy(MeV) 


Total(MeV) 


Exp(MeV) 


p 


1.52 


900 


201.34 


-141.46 


-67.21 


73.67 


966.34 


938 


A+ 


1.52 


900 


218.08 


-134.74 


143.01 


78.73 


1205.09 


1232 


A 


1.72 


900 


231.72 


-124.18 


99.50 


114.26 


1221.31 


1232 


A 


1.33 


1121.74 


185.36 


-147.26 


-77.16 


52.04 


1134.71 


1115 




1.44 


1121.74 


216.90 


-145.18 


-40.49 


61.41 


1214.37 


1194 




1.33 


1121.74 


185.36 


-147.26 


-42.80 


52.04 


1169.07 


1192 




1.44 


1121.74 


232.67 


-144.13 


139.96 


66.09 


1416.31 


1385 




1.33 


1121.74 


185.36 


-147.26 


179.94 


52.04 


1391.81 


1385 





1.51 


1565.22 


175.71 


-144.32 


28.75 


76.87 


1702.22 


1672 



Table 3: The full set of octet and decuplet particle energies in the TF quark model. 



VI. THREE CASE STUDIES 

In this section we will apply the fit found in the previous section to begin to explore some of the 
interesting phenomenology of high number multi-quark states. In particular, we will examine three 
possibilities: the i/-dibaryon, nucleon-nucleon 6 quark resonances and high multi-quark strange 
states. 

A. //-Dibaryon Considerations 

The i!f-dibaryon, which is an isospin / = 0, total angular momentum J = uuddss flavor state, 
is an interesting application of the model. Although there are many lattice results [THE], there are 
systematic problems in the calculations due to lattice volume and quark mass effects. There are 
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experimental results [151 [16] which limit the bomid state energy to shallow values if indeed it exists 
at all[17J. The lattice results must be extrapolated in quark mass, and the results of this procedure 
are still rather uncertain[7j. Our model is designed to evaluate states for further lattice study, so 
it is interesting to study this possible bound state with the new tool we now possess, even though 
we are doing this study "post-lattice" at this point. 

We may form a total I = 0, J = combination from combining two particles in the (AA), (SS) 
or (EN) systems[5j. The individual particle isospins are Ia = 0, /s = 1, /at = ^ and /= = 5- 
One can look up the Clebsch- Gordon coefficients to couple these systems to a \Itot = > state. 
In quantum mechanics the states mix with one another, giving a coupled channel problem[H]. 
Such considerations are beyond the present model's purview. However, given the quantum wave 
functions each of these combinations may be examined separately for the lowest mass. 

The possible TF quark 4 quark uuddss J = combinations are only three in number (we don't 
list some redundant configurations in the m„ = isospin limit): 

{u^u^d^d^s^s^) : Ni = 2,gi = 2,= 2;Ns = l,gs = 2, 
(u^uUUhh^) ■.Ni = 2,gi = 1;N2 = 1,92 = 2;Ns = 2,g, = 1, 
(u^u^dUh^s^) : iVi = 1,51 = 4;Ns = 1,5. = 2. 

Unfortunately, the u,d asymmetric second combination, {u^u^d^d'^s'^ s"^), can not be formed here 
because it involves three separate and different quark wave functions. This limits us to investigation 
of the AA state only, although this is also likely the lightest state from Pauli blocking considerations. 

Our result for the A A / = 0, J = state is 2360 MeV/c^. The component energies and radius 
are: 

Quark rest mass = 2244MeV/c^, 
Kinetic energy = 356MeV/c^, 
Potential energy = — 266MeV/c^, 
Volume energy = 85.3MeV/c^, 
Spinenergy = — 58.5MeV/c^, 
Radius = 1.62 fm 

This is 90 MeV/c^ more than the AA threshold. Our present results indicate the i7-dibaryon state 
is not bound in the TF quark model. 
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B. Nucleon-Nucleon Resonances 

When bound states are not present in a channel, our model is useful in predicting multi- 
quark baryon resonances. In Fig. 1 we show a "scan" of the configuration energies for the two 
wave function case for nucleon-nucleon scattering. Note these are the configuration energies, not 
the particle energies. To obtain particle resonance energies, we must proceed by projecting the 
quantum state into the particle configurations. Linear combinations of these configuration energies 
will be the actual resonance energies. 

It turns out there are 8 different configuration in this case. Moving from left to right in the 
figure displays how the energy rises as more particle are put into a given flavor state. Given that 
+ 52-^2 = 6, the explicit cases are: 





= 0,;iV2 


= 2,92 = 


= 3; 'Wi.iV2" = 0.2 




Ni 


= 0,;iV2 


= 3, 92 = 


= 2; = 0.3 




iVi 


= 0,;iV2 


= 6, 92 = 


= 1; "ATi.ATa" = 0.6 




iVi 


= 2,91 = 


2;A^2 = 


1,92 = 2; "iVi.Ar2" 


= 2.1 


iVi 


= 3,91 = 


i;iV2 = 


1,92 = 3; 'Wi.iV2" 


= 3.1 


iVi 


= 4,91 = 


i;A^2 = 


1,92 = 2; 'Wi.iV2" 


= 4.1 


ATi 


= 4,91 = 


i;A^2 = 


2,92 = l;"iVi.iV2" 


= 4.2 


iVi 


= 5,91 = 


l;iV2 = 


l,92 = l;'Wi.Ar2" 


= 5.1 



Since each configuration is uniquely specified by the iVi and N2 values, we use these to numbers to 
form a fractional number defined by N1.N2 to plot the results in Fig. 1. The integer number rep- 
resents the number of A^i particles and the decimal number represents the number of N2 particles. 
For example, 0.6 on the N1.N2 axis means Ni = 0; N2 = 6, g2 = Note that we only calculate the 
bare mass here in this quick scan of configuration energies. 

The lowest configuration energy is seen to be in the N1.N2 = 0.2 state at 2016 MeV/c^. This 
is as we expect from the exclusion principle since we are spreading 6 quarks across the maximum 
number of light flavors. In addition, the highest mass is N1.N2 = 0.6 at 2294 MeV/c^. This is 
also as expected. We see that all of the states are greater than the two nuclcon threshold and that 
there are no elemental six quark nucleon-nucleon bound states in our model. There is a resonance 
above two nucleon threshold in both pp and np scattering, but is impossible to tell if the states 
contributing are mesons and baryons or true 6 quark states. 
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FIG. 1: Nucleon-nucleon six quark configuration resonance chart. See the text for the meaning of the 
"iVi.A^2" notation. 

C. High Multi-Quark Strange States 



We have studied the possibiUty of a bound if-dibaryon in Section VIA Lattice results have 
already begun to accumulate for this system. We can get ahead of the lattice results by beginning 
to investigate the possibility of bound states with many more quarks. 

We have investigated the possibility of addition / = states constructed from A, A wave 
functions. Such states can be thought of as multi-combinations of ^/-dibaryons with 12, 18, 24 . . . 
quarks. Such states would be very difficult if not impossible for the lattice to simulate. Here, it is 
a simple matter of changing the initial simulation parameters to accommodate these new states. 
In this model, we formally have a maximum of 4 light quark flavors. We have investigated the 12 
quark configuration, 

Ni = 2,gi=A;Ns=2,gs = 2, 

the 18 quark configuration. 



Ni = 3,gi = A;Ns = 3,gs = 2, 
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and the 24 quark configuration, 

Ni=4,gi=4;N, = 4,gs = 2. 

Higher quark number states are obvious generalizations. We find the fohowing energies and masses 
in our investigation: 

For the 12 quark case, we find the mass as 4777 MeV/c^, 237 MeV/c^, more than the 4 A threshold 
and 57 MeV/c^ more than twice the TF ^f-dibaryon mass. The component energies and radius 
are: 

Quark rest mass = 4487MeV/c^, 
Kinetic energy = 596MeV/c^, 
Potential energy = — 443MeV/c^, 
Volume energy = 159MeV/c^, 
Spin energy = -22.0MeV/c^, 
Radius = 1.92 fm 

For the 18 quark case, we find the mass as 7228 MeV/c^, 417 MeV/c^ more than the 6 A threshold 
and 147 MeV/c^ more than three times the TF i7-dibaryon mass. The component energies and 
radius are: 

Quark rest mass = 6731MeV/c^, 
Kinetic energy = 861MeV/c^, 
Potential energy = — 601 MeV/c^, 
Volume energy = 251MeV/c^, 
Spin energy = -13.1 MeV/c^, 
Radius = 2.24 fm 

Finally, for the 24 quark case, we find the mass as 9679 MeV/c^, 591 MeV/c^ more than the 8 A 
threshold and 239 MeV/c^ more than four times the TF ^/-dibaryon mass. The component energies 
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and radius are: 

Quark rest mass = 9874MeV/c^, 
Kinetic energy = 1122MeV/c^, 
Potential energy = -751MeV/c^, 
Volume energy = 344MeV/c^, 
Spin energy = — 9.4MeV/c^, 
Radius = 2.49 fm 

These states are moving further away from the bound particle thresholds. Interestingly, the systems 
are remaining relatively non-relativistic, with a smaller percentage of the component energy being 
kinetic as more quarks are added. In addition, the spin energies are becoming a smaller system 
component as well. 

VII. CONCLUSIONS 

We have begun to develop the TF quark model as a tool in strong interaction dynamics. We 
have shown how to introduce spin into the TF quark model and found parameters which effectively 
fit the low energy baryon states. The results presented here show how to construct the necessary 
wave functions in the two wave function multi-flavor case with quantum numbers Ni, gi, N2, g2- 

We have undertaken three case studies: the ii-dibaryon, nucleon-nucleon resonances and high 
multi-quark strange states. Our results for strange quark states show growing systems which are 
remaining relatively non-relativistic and whose spin interactions are a decreasing percentage of the 
total energy. There are no hints that the systems can be bound. However, these results should be 
regarded as initial indictions rather than final results. Wc have also confirmed a phenomenology 
in nucleon-nucleon scattering consistent with known low energy hadronic physics. 

There are many remaining areas in which to continue the development. We note that the 
model may be extended: 

• To be relativistic 

• To include anti-quarks 

• To include central heavy quarks 

• To look at exotic forms of matter (eg., "strange" matter, color-flavor locked quarks from Cooper 
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pairs, or massive gluon models, or combinations) 

The relativistic equation to be solved is given in Ref . [Tj , although the system energy formulas are 
not yet developed. We can not strongly defend our non-relativistic assumption other than pointing 
out that the present model should be regarded as a work in progress. The inclusion of anti- 
quarks would simply be accomplished by introducing a different color coupling for such particles 
and averaging over the couplings between quarks and anti-quarks. Heavy central quarks can be 
examined using Coulombic sources much like atomic TF systems. Perhaps most interestingly, the 
color-flavor locking scenario, partly examined in Ref.jTj can also be examined. 

We believe that the TF quark model, in combination with lattice studies, can be an effective 
"tool" in evaluating characteristics of interesting baryonic states. 
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